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We investigate efficient entanglement purification through double selection process. This method 
works with higher noise thresholds for the communication channels and local operations, and achieves 
higher fidelity of purified states. Furthermore it provides a yield comparable to the usual protocol 
with single selection. It is shown by general considerations that the double selection is optimal to 
remove the first-order errors, achieving the upper bound on the fidelity of purified states in the low 
noise regime. The double selection is also applied to purification of multi-partite entanglement such 
as two-colorable graph states. 



PACS numbers: 03.67.Hk, 03.67.Pp, 03.67.-a 



I. INTRODUCTION 

Recently a number of protocols based on entangle- 
ment have been developed in quantum communication 
and computation. For example, quantum teleportation, 
superdense coding, quantum cryptography and quantum 
repeater employ bipartite entanglement [l|, 0, H, 0] ■ Clus- 
ter state computation, quantum error correction and 
multiparty cryptography utilize multi-partite entangle- 
ment 0, Si 0| ■ "^^^ performance of these entanglement- 
based protocols highly depends on the fidelity of entan- 
gled states. That is, high fidelity entangled states are 
essential for secure communication and reliable compu- 
tation. In this viewpoint, it is a very important task to 
prepare and share high fidelity entangled states. 

A way to share high fidelity entangled states via noisy 
communication channels is known as entanglement pu- 
rification [1, . This protocol is proposed originally to 
share EPR states, and then extended for a large class of 
multi-partite entangled states, including the GHZ states, 
two-colorable graph states, stabilizer states and W states 
[13, [m, [l^ Il4l. [l5l | . In a bipartite situation with noisy 
channels and perfect local operations, we pre-purify ini- 
tial states with a recurrence protocol, which gives a high 
threshold for the noise of the communication channel but 
a low yield of purified states. Then, a hashing protocol 
may be implemented to get pure entangled states with 
a nonzero yield. The hashing protocol, however, breaks 
down as soon as local operations become slightly imper- 
fect. The entanglement purification in such a situation 
is first analyzed in the context of quantum repeater [3] , 
where the usual recurrence protocol [1, [§] is applied with 
imperfect local operations. In this situation, the fidelity 
of purified states is limited, and noise thresholds are re- 
quired for the local operations as well as the communi- 
cation channels. This is the most distinct point from the 
case with perfect local operations. Thus, in order to real- 
ize entanglement-based protocols using practical devices, 
which inevitably have imperfections, we need to develop 
efficient purification methods, which work well even with 
noisy communication channels and imperfect local oper- 
ations. 

In this paper we investigate an efficient entanglement 



purification protocol with more accurate post-selection 
through double verification process. This method with 
double selection enables us to achieve considerably high 
fidelity of purified states even with actual imperfection. 
Furthermore it works with higher noise thresholds for the 
communication channels and local operations compared 
with the single-selection scheme 0, |9|. It may be con- 
sidered that the elaborate post-selection consumes more 
resources to decrease the yield of purification. However, 
this is not necessarily the case. The present method 
with double selection provides a yield comparable to the 
single-selection scheme, especially when higher fidelity 
states are required with noisier local operations. It is 
really shown that the double selection is optimal to re- 
move the first-order errors, achieving the upper bound on 
the fidelity of purified states in the low noise regime. The 
double selection is also applied to purification of multi- 
partite entanglement such as two-colorable graph states. 
As an example, we numerically investigate the perfor- 
mance for a CSS code (Steane's 7-qubit code), and com- 
pare it with a multi-partite purification protocol with 
single selection [l3|. As the double selection can achieve 
considerably high fidelity under noisy channels and oper- 
ations, it is indeed profitable for quantum computation 
as well as quantum communication. 



The rest of the paper is organized as follow. In Sec. 
|TT]we investigate the double selection in the bipartite en- 
tanglement purification. The performance of the double- 
selection protocol is analyzed in detail and compared 
with the usual protocol with single selection. The op- 
timality of the double selection is also discussed in the 
low noise regime. The double selection is then applied to 
the multi-partite entanglement purification in Sec. IIIIl 
where as an example the purification of the Steane's 7- 
qubit code is numerically investigated. Sec. llVl is devoted 
to conclusion. 
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FIG. 1; Bipartite entanglement purification with single selec- 
tion. 



II. BIPARTITE ENTANGLEMENT 
PURIFICATION 

A. Single selection 

We first review the usual recurrence protocol for pu- 
rification with single selection [1, |^ . Suppose Alice and 
Bob share EPR pairs, which are decohered via a noisy 
quantum channel, and purify them by using noisy C-Not 
gates and a classical channel. The recurrence protocol 
[1, consists of two copies of an EPR pair, a bilateral 
C-Not gate and a bilateral measurement at each round 
of purification (see Fig. [1]). Here, a bilateral operation 
means a tensor product of two identical operations which 
are simultaneously implemented by Alice and Bob. The 
purification procedure is specifically described as follows: 

(i) Alice and Bob share two identical EPR pairs p'*'-' 
and p^^^ through a noisy quantum channel. 

(ii) They operate a bilateral C-Not gate for p^"^ and p'-'^' 
as the control and target qubits, respectively. 

(iii) They bilaterally measure p'^' in the Z basis 
{|0),|1)}, and obtain the measurement outcomes 
nia (Alice) and rrih (Bob). 

(iv) They communicate these measurement outcomes to 
each other. Then, they keep p^'^^ if the measurement 
outcomes coincide as TOq = 'Tif,. Otherwise, they 
discard p'"'. 

Alice and Bob iterate the procedures (ii)-(iv) by using the 
output states which survive the selection in (iv) as the in- 
put states for the next round of purification where the X 
and Z bases of their reference frames are exchanged with 
a Hadamard transformation i. As seen in the above, a 
single bilateral operation determines whether p^^^ should 
be kept or discarded at each round. That is, the purifi- 
cation procedure is made with single selection. 

The noisy EPR pairs p'-"-' and p^^^ are given as two 
copies of a Bell diagonal state p, 
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(cto = I)- In the case with ideal (perfect) local opera- 
tions, the above purification procedure with single selec- 
tion provides the output states of p^"^ depending on the 
initial states of p'^' as 

X 

where x denotes the output states discarded by the post- 
selection. Alice and Bob are to distill \ 4>q) with the largest 
probability Fq (with a suitable basis transformation if 

necessary). For the most likely input state (^q""* (g 4>\^'^ 
with the probability -Fq , the purification procedure gives 
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the desired output state 0o ■ other hand, the 

input state cfyf^ ® <j)"^^ with the probability FqF^ results 
in a wrong output state (j)"^^ , and so on. This may be 
viewed as the 0-3 error of (j)'^'' = (0-3 ® o-o)0o^^(o'3 ® ctq) 
in p(^^ is propagated to p*-°^ through the purification pro- 
cedure. In this way, the purification procedure provides 
a map p ^ p' for the Bell diagonal states. This is even 
the case when the Pauli noise is introduced for the lo- 
cal operations, which is described later. This purifica- 
tion map can be represented by a i?^ map S from 
the input state- vector F — (Fq, ^1,^2,^3) to the output 
state- vector F' = (F(5, Fj', F^, F^) as 

F'-5(F), (4) 
or in terms of a transition probability tensor Sf^ 



as 



F' = 



(5) 



where 



P5(F)=^SfF,F, (6) 

i'jk 

is the success probability responsible for the normaliza- 
tion "^^Fl — 1. The transition probability tensor is 

i 

given including the noise parameters for the local oper- 
ations and the exchange of the reference frames at each 
round (see Appendix A for the detail). Then the maxi- 
mum achievable fidelity of purified states is determined 
as an appropriate fixed point of the map S. 

B. Double selection 

We now present the entanglement purification with 
double selection. The double-selection protocol consists 



3 



Alice 




^ ( 


1- 




^ € 


H 


H 



Bob 



FIG. 2: Bipartite entanglement purification with double se- 
lection. 



of three copies of EPR pairs, two bilateral C-Not gates 
and two bilateral measurements at each round of purifi- 
cation (see Fig. [2]), as described in the following: 

(i) Alice and Bob share three identical EPR pairs p^"^ , 
/O^^^ and p^^-* through a noisy quantum channel. 

(ii) They operate a bilateral C-Not gate for p^*^^ and p^^^ 
as the control and target qubits, respectively. 

(iii) Next they operate a bilateral C-Not gate for p*^^^ and 
p^^^ as the control and target qubits, respectively. 

(iv) They bilaterally measure p^^^ and p*^^^ in the Z and 
X bases respectively, and obtain the measurement 



outcomes mi^\ mi^'' 



(Alice) and m 



(1) 

b 



m 



(2) 



(Bob). 



(v) They communicate these measurement outcomes to 
each other. Then, they keep p'-^-' if both of the mea- 



surement outcomes coincide as nia 



(1) 



,(2) 



,(2) 



m'^i^^ and 



Otherwise, they discard p^"-* . 



Similarly to the single-selection protocol, Alice and Bob 
iterate the procedures (ii)-(v) by using the output states 
which survive the selection in (v) as the input states for 
the next round where the X and Z bases of their reference 
frames are exchanged with a Hadamard transformation. 

The action of these procedures with ideal local opera- 
tions is described as 
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with p(^) — and so on with p^^^ — 

The double selection removes the errors more efficiently 
by using two ancilla EPR pairs. For example, as seen in 
the above table for (/)q°'' (g) (f)^^"^ with ^q^-* , we can detect 
and discard the Z error propagation as 0q°'' (g) (f)'^^ (/ig"' 
which is the first-order failure event in the single-selection 
scheme. Thus, this double selection process significantly 



improves the maximum achievable fidelity of purified 
states. The optimality of the double selection to detect 
the first-order errors will be discussed in detail later. 

The purification map V of the double selection is de- 
scried as 



F' - P(F), (7) 
or in terms of a transition probability tensor Dj''^ as 



jki 



where 



pr,iF) = Y,Dt'F,F,F^ 

ijkl 



(9) 



is the success probability responsible for the normaliza- 
tion. The transition probability tensor D^*^' is given in 
Appendix A. 



C. Performance analysis and comparison 

For each scheme, we investigate the performance by 
considering the working range of the local operations, the 
maximum achievable fidelity of purified states, the min- 
imum fidelity required for the quantum communication 
channels and the EPR resources consumed to achieve a 
target fidelity. Here we assume that the imperfect C-Not 
gate is implemented as a perfect C-Not gate followed by 
two-qubit deporalizing errors with probabilities pij : 



Af{p) = {l-pg)p+ ^ pij{ai®(jj)p{(7i®(7j), 
(»j)/(o,o) 



(10) 



where i,j = 0, 1,2,3, and pg — pij. The imper- 

(»j)/(o,o) 

feet measurement in the Z basis is described with POVM 
elements and error probability Pm as 

^0 = (1-p™)|0)(0H-p™|1)(1|, (11) 
El = (l-p„0|l)(l|+p™|0)(0|. (12) 

The imperfect measurement in the X basis is described 
similarly as 



HEqH, 
HEiH, 



(13) 
(14) 



where H denotes a Hadamard operation. Then the tran- 
sition probability tensors Sf'' and Df'^'' , which charac- 
terize the purification maps S and V respectively, are 
calculated including the error probabilities pij and Pm 
(see Appendix A). By operating the purification map A 
recursively, we can calculate the state-vector F^"^ at the 
n round of purification iteration as 



(15) 



4 




E 
Q. 




0.1 0.12 



FIG. 3: Purification curve. The fidelity of the output state 
F' is plotted as a function of the fidelity of the input state F. 
Curves from top to bottom correspond to perfect local oper- 
ations, imperfect local operations inside the working range, 
and those outside the working range. 



FIG. 4: Working range {pg,Pm) of local operations. Each 
purification map A has the non-trivial fixed points Fmax and 
-FiTiin for the error probabilities (pg,Pm) below the threshold 
curve (solid line for T) and dotted line for S). 



{A = S for the single selection and T) for the double se- 
lection) . The EPR states are shared initially through the 
noisy communication channel C as 

^ ^ 

C{p)^Fchp-\ (To pat Q-Q. (16) 

^ ^=l 

Thus we set the initial state- vector for C{(j)o) as 
1 — Fch 1 — ^ch 1 — ^ch 



F(o) = F, 



ch 1 



(17) 



In the case with ideal local operations, the purification 
map A has three fixed points for F = Fq and F' = Fq, 
i.e., the maximum achievable fidelity i^max = 1, the min- 
imum required channel fidelity -Fmin = 1/2 and the com- 
pletely mixed state i^mix — 1/4, as seen from the purifi- 
cation curve in Fig. [3l Then the imperfection of the 
local operations shifts down the purification curve. If 
the error probabilities Pg and pm are inside the work- 
ing range of the purification scheme, there still exist the 
three fixed points Fmax < 1, Fmin > 1/2 and Fmix = 1/4. 
In this case, if the communication channel has a fidelity 
Fch higher than Fmin, we can achieve the fidelity Fmax 
by iterating the purification procedure. On the other 
hand, outside the working range, the purification pro- 
cedure converts the noisy EPR state to the completely 
mixed state (Fmix = 1/4) with no entanglement. 

The working range {pg,Pm) of the local operations is 
shown in Fig. 0] for each purification scheme. Here the 
error probabilities of the C-Not gate are taken equally 
as Pij = Pg/15. As seen in Fig. [H the double-selection 
scheme has higher thresholds for the error probabilities 
than the single-selection scheme. Thus it works well with 
noisier local operations. If we choose a different error 
distribution pio = poi = g^, Pij = QiQj {i ^ j and i,j^O) 



T3 



1 

0.95 

0.9 
0.85 

0.8 
0.75 

0.7 
0.65 

0.6 
0.55 

0.5 

















\ 

\ - 




'"mirr/'^ 




' Double 


1 1 


Single 

1 1 ^1 1 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 
P=Pq=Pm 



FIG. 5: Maximum achievable fidelity Fmax (upper curves) 
and minimum required channel fidelity JVnin (lower curves) 
are plotted as functions of the error probability p = pg — pm 
for the double and single selections. 



and Pg 
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Qi, as adopted in Ref. [3], we estimate the 



threshold values 3 — 4% and 4 — 5% of Pg with pm = 
for the single and double selections, respectively, which 
are consistent with an upper bound 5.3% derived in Ref. 

m- 

In Fig. O Fmax (upper curves) and Fmin (lower curves) 
are plotted as functions of the error probability p {p — 
Pg — Pm for definiteness). The double selection achieves 
a higher fidelity of purification Fmax with a less fidelity of 
communication channel Fmin, compared with the single 
selection. These improvements are due to the fact that 
the double selection can detect more efficiently the first- 
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FIG. 6: The yield Ya{F, Fch = 0.8) is plotted as a function of 
the target fidelity F for each protocol with pg = pm = 0.02 
(upper curves) and Pg = pm = 0.04 (lower curves). 



order errors, as discussed in detail later. 

We compare the yields of the purification protocols. 
The yield Ya{F, Fch) is defined as the inverse of the num- 
ber of EPR pairs consumed to achieve a target fidelity 
F under the channel fidelity Fch It is calculated for 
each scheme A — S.V as 



YA{F,Fch) = 



(18) 



where «^(f, fch) denotes the minimum number of 
rounds which is required to achieve the fidelity F, 
p^(F^"~^'') denotes the probability to pass the purifica- 
tion procedure at the round n, as given in Eqs. ^ and 
([9]) , and Nj, denotes the number of EPR pairs consumed 
at each round {Ng — 2 and iVp = 3). We plot in Fig. 
[6] the yield Yj\{F,Fch — 0.8) as a function of the target 
fidelity F for each protocol with pg — = 0.02 (upper 
curves) and Pg — Pm = 0.04 (lower curves). With less 
noisy local operations {pg — pm — 0.02), both schemes 
provide comparable yields to achieve a target fidelity 
F w 0.9, where the numbers of rounds are = 4 (sin- 
gle) and n-p = 2 (double). On the other hand, even when 
noisier local operations {pg — — 0.04) are used, the 
double-selection scheme still provides a reasonable yield 
to achieve a target fidelity F « 0.9, where = 16 and 
n-p = 4. Since the double-selection scheme uses three 
EPR pairs at each round, it may be thought to cost 
more resources than the single-selection scheme with two 
EPR pairs at each round. However, this is not the case. 
As shown in the above, the double-selection scheme pro- 
vides a comparable or even better yield. This is because 
the double-selection scheme increases the fidelity of EPR 
pairs considerably faster than the single-selection scheme. 
In the case of the double selection, an additional EPR 
pair is worth enough to enhance the fidelity of purified 
states, as discussed in the following. 
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FIG. 7: Setup of a purification protocol. An upper bound on 
the fidelity is determined in the first order by the unavoidable 
errors (indicated by black stars) introduced by the final two 
C-Not gates. Similar bounds are obtained with other config- 
urations of two-qubit gates. 



D. Optimality of the double selection 

We here discuss the optimality of the double selection 
in the first order of the errors. We first consider an upper 
bound on the fidelity of purified states by general argu- 
ments. Then, it will be shown that the double selection 
really saturates this upper bound in the low noise regime. 

While purification protocols contain many imperfect 
local operations, we can reasonably derive an upper 
bound on the fidelity in the first order by noting the er- 
rors introduced by the final two gate operations applied 
on the source EPR pair. Here we assume that the local 
operations for purification are implemented by bilateral 
two-qubit Clifford gates and measurements. (One-qubit 
gates may be absorbed into relevant two-qubit gates and 
local reference frames.) Specifically, consider a purifica- 
tion protocol with a combination of two bilateral C-Not 
gates at the final step, as shown in Fig. [71 We inspect 
these final C-Not gate operations in the Alice's site to ob- 
serve the undetectable errors left on the output state of 
the source qubit (s). (The same argument is made in the 
Bob's site.) Passing through this set of the final gates, the 
preceding errors on the source qubit (s) are propagated 
to either or both of the ancilla qubits (al) and (a2). Thus 
these preceding errors are all regarded to be detectable, 
that is an optimal purification protocol is expected to re- 
move them by post-selection after the measurements of 
the ancillas. The upper bound on the fidelity is rather 
determined by some of the errors introduced by the final 
two gates themselves (black stars in Fig. [7]) which are 
undetectable and thus unavoidable without leaving any 
information on the ancillas. As for the second final C- 



Not gate, only the a^'^ (E) crg"^^ error with the probability 
P30 is undetectable, since it commutes with the final C- 
Not gate. On the other hand, all the (8) Tq errors 
of the final C-Not gate with the probabilities Pio are un- 
detectable inevitably, since the output source qubit does 
not interact with ancillas afterward. By adding these un- 
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detectable errors due to the final two C-Not gates in the 
sites of Alice and Bob, the upper bound on the fidelity is 
placed in the first order of the errors as 

Fuppcr = 1-2 |^P30 + J2p^0^ - ^^pI^ (19) 

Note here that the measurement error is not involved 
in considering the detectability of the errors with their 
propagation to the ancilla qubits. A portion of the right 
output state 4>q may be discarded by the errors in measur- 
ing the ancillas for verification. This slight reduction of 
the right state is cancelled by the rcnormalization after 
the post-selection except for the higher-order contribu- 
tions. Thus the measurement error does not contribute 
to the upper bound in the first order. The precise upper 
bound will be lower than this first-order bound due to 
the higher-order contributions of the gate and measure- 
ment errors, though it is beyond our scope to estimate 
them analytically. 

Another protocol may be considered by exchanging the 
final two C-Not gates in Fig. [71 (This is actually the case 
in the recurrence protocol considered so far when the pu- 
rification procedure is finished at an even round.) Simi- 
larly, by observing the undetectable errors introduced by 
the final two C-Not gates, the upper bound on the fidelity 
is placed as 

^uppcr = 1 - 2 l^pio + Y.P^O^ - 0(Pg). (20) 

We note that these upper bounds in Eqs. ^9]) and ([20]) 
are valid even if the final two C-Not gates are set in the 
same direction (e.g., both the C-Not gates are controlled 
by the source qubit). In such configurations, some of 
the preceding errors on the source qubit become unde- 
tectable, commuting with the final two gates. Thus these 
configurations merely lower the fidelity. 

These two upper bounds i^upper and -Fuppcr coincide 
with each other for the uniform distribution of the gate 
errors pij = Pg/15, as adopted in calculating the maxi- 
mum achievable fidelities of the single and double selec- 
tions (Fig. [5]). In general, the upper bound is given by 
max[Fiippor, ^uppcr] depending on the error distribution; 
in a recurrence protocol we should determine whether 
the purification procedure is finished at an even or odd 
round. We may also use other two-qubit Clifford gates 
instead of C-Not gates. Then, we obtain a similar up- 
per bound with a suitable permutation among pio's in 
Eq. or Eq. (|20p . by considering which error of the 
second final gate is left commuting with the final gate. 

We now show that the double selection really removes 
all the detectable errors in the first order, achieving the 
upper bound F^pper- (We here adopt the uniform error 
distribution for definiteness.) The recurrence protocol 
considered so far, with either single or double selection, 
has the setup as shown in Fig. [7] with the proper ex- 
change of the reference frames at each round. In the 



single selection (Fig. [T]), the cti and (T2 (o'2 and 0-3) er- 
rors on each ancilla EPR pair are detected by the Z (X) 
measurement, while the 173 (cti) error is not detected. In 
the double selection (Fig. the primary ancilla EPR 
pair is further verified with the secondary ancilla EPR 
pair so that all the errors on the primary ancilla EPR 
pair are detected in the first order. That is, all the de- 
tectable errors on the source EPR pair are certainly re- 
moved by the double selection through the interaction 
with the primary ancilla EPR pairs, as expected in de- 
riving the upper bound. Therefore, the double-selection 
protocol achieves the upper bound in the first order as 

FZ. = F„~0{pI). (21) 

This has been really confirmed by the numerical calcu- 
lation in the low gate noise regime Pg < 0.02, almost 
independently of the measurement error p„i < 0.05. 

On the other hand, as noted above, in the single se- 
lection (Fig. [T|) the CT3 (cti) error on the ancilla EPR 
pair cannot be detected by the Z (X) measurement at 
an odd (even) roimd. Even in this case, all the errors be- 
fore the second final C-Not gate are detected in the first 
order through the final two rounds, because the reference 
frames are properly exchanged at each round. However, 
some of the errors which have been excluded as detectable 
in deriving -Fuppor are actually left in the single-selection 
protocol. By counting these extra errors due to the final 
two C-Not gates in Fig. [3 we find that the maximum 
achievable fidelity of the single selection decreases as 

FL. = FZ. - (8/15)p, - 0{pI). (22) 

As seen so far, the double selection is necessary and 
sufficient to remove fully the detectable first-order er- 
rors, saturating the upper bound on the fidelity in the low 
noise regime. There is no room for the triple (or more) se- 
lection to improve the fidelity except for the higher-order 
error contributions. In practice, as the triple selection 
itself introduces more errors with additional noisy oper- 
ations, it will hardly improve the maximum achievable 
fidelity of the double selection. 

III. MULTI-PARTITE ENTANGLEMENT 
PURIFICATION 

Recently entanglement purification is applied to a 
large class of multi-partite entan glem ents including two- 
colorable graph states [ll|, [Tsl. H^. We can extend 
the double-selection scheme to such protocols for multi- 
partite entanglement purification. Specifically, we here 
consider the purification of two-colorable graph states 

m- 

A graph state is expressed as an eigenstate of a set of 
stabilizer operators 

Kj = X, (g) Zk, (23) 



7 



where Xj and denote the PauH operators acting on 
the j-th and fc-th qubits respectively, and Vj denotes the 
set of vertices which are connected to the j-th qubit [l3| ■ 
Then a graph state is described as 



KM,, 



(24) 



where (— = 0,1) denotes the eigenvalue ±1 of 
the operator Kj. Especially a graph state which can 
be divided into two partitions A and B such that no 
vertices within one set are connected by edges is called a 
two-colorable graph state. It is described as 



I t i t I 
t 1 t i t 
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FIG. 8; Purification of two-colorable graph states with double 
selection. 



(25) 

where /x^ and iig denote the eigenvalues of the vertices 
within the partitions A and B, respectively. 

The entanglement purification protocol with double 
selection for a noisy mixture of n-partite two-colorable 
graph states 

E v.,/^j*/^.,/^J(*/^.^mJ (26) 

is implemented as follows (see Fig. |8]): 

(i) Alice, Bob, • • • , Nancy share three identical two- 
colorable graph states p^^\ p^^^ and p^^^ through a 
noisy quantum channel. 

(ii) They operate a multi-lateral C-Not gate for p^^"^ and 

with the partition A [B] of [p^^)] as the 
control and the partition A [B] of p'^^ [p''^^] as the 
target, respectively. 

(iii) Similarly they operate a multi-lateral C-Not gate 
for p(2) and p^^\ 

(iv) They multi-laterally measure the partition A [B] of 

and the partition A [B] of p^^-* in the Z and 
X [X and Z] bases respectively, and obtain mea- 
surement outcomes (—1)^'* ' [("1)"°^' '] and (—1)^'* ' 

[(—1)^3 ], where £,i,Cj ~ 0,1, and the qubit i [j] 
belongs to the partition A [B] . 

(v) They communicate these measurement outcomes to 
each other. Then they keep p'-^-' if for all i and 

J, Ci'^ ® E.ey, ei'^ = and e Efcey. ^ = 0, 
which implies /x^^®/.i^^©/.i^^ = and /x^^©/x^^ = 
0, respectively, where ® denotes bitwise addition 
modulo 2. 

They iterate the procedures (ii)-(v) by using the out- 
put states which survive the selection in (v) as the input 
states for the next round where the X and Z bases of 
their reference frames are exchanged with a Hadamard 
transformation, similarly to the bipartite case. These 



procedures with perfect operations act as 



^,,(0) ,,(0)/ 






(27) 








(28) 








(29) 



Thus in the double-selection scheme we can detect even 
the propagation of the errors /it^'' from p'^^^ to p^^'' by 
using p^2) the same as the bipartite case. 

As an example of the multi-partite purification with 
the double selection, we numerically investigate a CSS 
code state, specifically the Steane's 7-qubit code state, 
which is a two-colorable graph state. Then i ts p erfor- 
mance is compared with the ADB protocol [11|. We 
consider the multiparty communication situation, where 
n-qubit two-colorable graph states l^o^.Ois) (^Oa.Ois | are 
shared through n identical noisy channels C®". The noisy 
copies Pin = C'^"{\^0a.0b){'^0a,0b\) of the n-qubit en- 
tangled state are purified with the noisy C-Not gates 
and measurements the same as the bipartite case. We 
simulate directly the noisy operations in the communi- 
cation and purification procedures by using Monte-Carlo 
method. (It is highly complicated to provide the purifi- 
cation map in terms of the transition probability tensor.) 
The fidelity of the purified state p' is measured by 

Fip', |*o.,oJ) = (*o,,o«|p'|*o.,oJ. (30) 
If the initial fidelity 

^ i^(pi„, I^Oa.Os)) = + 0{{1 - Feh)') (31) 

is higher than Fmin, we can achieve the fidelity F^ax by 
iterating the purification procedure. 

The resultant maximum achievable fidelity -Fmax and 
minimum required initial fidelity i^min are plotted for the 
Steane's 7-qubit code state 1^*0^, Ob) — |0i) as functions 
of the error probability p = Pg ^ Pm in Fig. [S] As we 
expect, the double-selection protocol achieves consider- 
ably higher fidelity of the purified states requiring lower 
fidelity of the channels. Then the noise threshold for the 
local operations is improved from 5—6% (single) to 8—9% 
(double). It is also seen in Fig. (TUlthat both schemes pro- 
vide comparable yields, similarly to the bipartite case. 
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profitable also in quantum computation. In fact, en- 
coded ancilla qubits are used to stabilize a computation 
in a fault-tolerant way, and the performance of computa- 
tion highly depends on the fidelity of these ancilla qubits 
[H, Hal- In the usual fault-tolerant context [H, Ho, HI, 
these encoded ancilla qubits are prepared through the 
single selection. Thus the double selection has a good 
potential to improve the noise threshold of fault-tolerant 
computation. The verification process with the double 
selection is used in fault-tolerant computation with con- 
catenated construction of verified logical cluster states, 
achieving a considerably high noise threshold ~ 4% [l^ . 
We will discuss elsewhere further applications of the dou- 
ble selection for quantum computation. 



FIG. 9: Maximum achievable fidelity i^max (upper curves) and 
minimum required diannel fidelity Fmin (lower curves) for the 
Steane's 7-qubit code state |0l) are plotted as functions of 
the error probability p = pg — pm for the single and double 
selections. 
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FIG. 10: The yield y^(_F,Fch = 0.9) is plotted as a function 
of the target fidelity F for each protocol with pg — pm ~ 0.02 
(upper curves) and pg = pm = 0.04 (lower curves). 



IV. CONCLUSION 

We have investigated entanglement purification with 
double selection. It has been shown that the maxi- 
mum achievable fidelity can be improved significantly by 
the double selection under noisy communication channels 
and imperfect local operations. Furthermore the double- 
selection scheme provides a reasonable yield, which 
is comparable or even better than the single-selection 
scheme. The double selection is really optimal to remove 
the first-order errors, achieving the upper bound on the 
fidelity of purified states in the lower noise regime. These 
results really indicate that the double selection is more 
suitable for realization of entanglement-based protocols. 
The double selection is also applied to multi-partite en- 
tanglement purification, specifically two-colorable graph 
states. The improvement of the fidelity has been shown 
numerically for the Steane's 7-qubit code state. The dou- 
ble selection can be extended to all graph states in the 
same way. Since multi-partite entangled states, such as 
CSS codes and cluster states, play very important roles 
in quantum computation as well as quantum communica- 
tion, the double selection has a good potential to improve 
the performance of quantum computation. 



The yields are, however, significantly low compared with 
the bipartite case. This is because the coincidence of the 
more measurement outcomes is required in the multi- 
partite case, reducing the success probability of post- 
selection. The upper bound on the fidelity with imperfect 
local operations can be extend to the two-colorable graph 
states in the same line. Then the maximum achievable 
fidelity Fmax of the double-selection scheme saturates the 
upper bound Fuppcr = 1 — 7 x (4/15)pg for n = 7 with 
Pij — pg/15, which is really confirmed in the low noise 
regime by the numerical simulation. 

The two-colorable graph states, including cluster and 
CSS code states, play very important roles in quantum 
computation as well as quantum communication. Then 
these results really indicate that the double selection is 
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APPENDIX A: DERIVATION OF TRANSITION 
PROBABILITY TENSORS 

The purification procedure of the single selection be- 
fore the post-selection is described as a linear map S of 
the two Bell states c/jf ^ and (f)'j^^ as 

5(01'^) «,/,«) = 0(1). (Al) 
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This map consists of the noisy bilateral C-Not gate The bilateral Z-basis measurement is given by 



surement with error M ( 
of the reference frames at each round by the bilateral H 
transformation (ideal for simplicity) W(<;&i°^) = H^cj)'^^: 



G 



the bilateral Z-basis mea- 
: Ml'(j)[^^ and the exchange 



ah' 



The noisy bilateral C-Not gate G^-'^ is decomposed into 

the ideal bilateral C-Not gate [/^■^ and the bilateral com- 
bination N^f^N^^ of the C-Not gate noises at Ahce and 
Bob: 



(A3) 



Note here that the Bell states are two-colorable graph 
states up to a local Hadamard operation as 



(A4) 



where 



= (Mil, m's) = (0, 0), (1, 0), (1, 1), (0, 1) (A5) 

for i = 0,1,2,3, respectively. In the graph state repre- 
sentation, the action of the ideal bilateral C-Not gate is 
described as 



1* 



(0) 



or simply 

ZY(/x* (g) fj,^) 



I* 



. , (A6) 



{ji\,y^B®i^B)® (Pa © > A* B ) • ( A7) 



Then the transition tensor of the ideal bilateral C-Not 
gate is given by 



ab 




-- Uin" (g fj.'')] 



(A8) 



e.g., Vll = 1 for /Lx2 = (l,l) (l,l) ^ (1,0) 
(0, 1) = fi^ ® (1^ by U, and so on. The noise on the 
C-Not gate is given by 

= P^l G V.^ , ibd) e v., ] , (A9) 

e.g., A^"" = pio for (ac) = (01) with ai®aQ(j)Qai®aQ = 
and (bd) = (00) with ctq ao^ofo ctq = ^o, and so on, 
where poo = 1 — Pij = ^ — Pg, and cTj's provide 

(iJ)/(0,0) 

the permutations among the Bell states ^j's as 



e™6 



= m^m: 



(A14) 



(A2) as a product of the single measurements given by 



m: = < prn [(be) G Va,] 

[{be)€V^„V^,] 



The bilateral Hadamard operation is given by 



(A15) 



^k = KK' 



(A16) 



where Hq 



hi 



hi 



hi = 1 and /i° = for the 
others. After all the transition probability tensor Sf'^ of 
the single selection is obtained by picking up the post- 



selected states with Z = 0, 3 from Sfl' as 



S: 



(A17) 



Similarly, the purification procedure of the double se- 
lection before the post-selection is described as a linear 
map V of the three Bell states (l)f^^ , 



bf\ 4^ as 



Imn 



(A18) 



where 



(A19) 




12 3 

10 3 2 

12 3 

2 3 1 



(AlO) 
(All) 
(A12) 



provides the X-basis measurement. Then the transition 
probability tensor of the double selection is obtained by 
the post-selection as 



I / J Iran 

m=0,3:n=0,l 



(A20) 



^(73 — 



12 3 
3 2 10 



(A13) 
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